New Superhard Carbon Phases Between Graphite and Diamond 



<N 

o 

<N 

a 

=3 



Chaoyu He, 1 L. Z. Sun, 1 ^] C. X. Zhang, 1 K. W. Zhang, 1 Xiangyang Peng, 1 and Jianxin Zhong 1 '!!) 

1 Institute for Quantum Engineering and Micro-Nano Energy Technology, Xiangtan University, Xiangtan ^11105, China 

(Dated: March 1, 2013) 

Two new carbon allotropes (H-carbon and S-carbon) are proposed, as possible candidates for 
the intermediate superhard phases between graphite and diamond obtained in the process of 
cold compressing graphite, based on the results of first-principles calculations. Both H-carbon 
and S-carbon are more stable than previously proposed M-carbon and W-carbon and their bulk 
modulus are comparable to that of diamond. H-carbon is an indirect-band-gap semiconductor with 
a gap of 4.459 eV and S-carbon is a direct-band-gap semiconductor with a gap of 4.343 eV. The 
transition pressure from cold compressing graphite is 10.08 GPa and 5.93 Gpa for H-carbon and 
S-carbon, respectively, which is in consistent with the recent experimental report. 

PACS numbers: 61.50.Ks, 61.66.Bi, 62.50. -p, 63.20. D- 
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INTRODUCTION 



Carbon is considered as the most active element in 
the periodic table due to its broad sp, sp 2 and sp 3 
hybridizing ability. Besides the four best-known carbon 
allotropes, graphite, cubic-diamond, hexagonal diamond 
and amorphous carbon, an unknown superhard phase 
of carbon has been reported in experiment jjj-lsl along 
with the structural phase transition in cold compressed 
graphite. Several structures have been proposed as 
the candidates for this superhard phase, such as the 
monoclinic M-carbon [6, 7], cubic body center C4 carbon 
(bct-C4) [8] and the orthorhombic W-carbon [9]. The 
cohesive energy of the most stable W-carbon is about 
160 meV per atom higher than that of diamond, indi- 
cating that other more stable carbon phases may exist. 
Very recently, another new carbon allotrope, named as 
Z-carbon, was proposed and investigated at almost the 
same time by three independent research groups [10l4l2] . 
Z-carbon is more stable (its cohesive energy is about 
129 meV per atom larger than that of diamond) and 
harder than W-carbon. Moreover, its transition pressure 
is around 10 Gpa which is lower than that of W-carbon. 
Thus, it is believed that the formation in the process of 
the cold compressing graphite prefers Z-carbon rather 
than W-carbon. Although none of them can fit the 
experimental results satisfactorily, these theoretically 
proposed phases are significant in understanding the 
experimental process of cold compressing of graphite. 

The formation of the new carbon phase needs large 
lattice distortions, such as buckling between adjacent 
carbon layers of graphite. Therefore, the phase tran- 
sition depends not only on the particular nucleation 
history [13| but also on the structure of the starting 
graphitic material. Intuitively, such buckling between 
carbon layers will start from the natural AB stacking 
graphite rather than the AA stacking graphite chosen in 
previous theoretical studies |7H12|. From this hypoth- 
esis, we propose two new superhard crystals, namely 




FIG. 1: Crystal structure of H-carbon(a), initial AB stacking 
graphite supercell for H-carbon (b) and side view containing 
five and seven carbon rings of H-carbon (c). Crystal structure 
of S-carbon (d), initial AB stacking graphite supercell for S- 
carbon (e) and side view containing five and seven carbon 
rings of S-carbon (f). 



H-carbon and S-carbon, as potential candidates for the 
experimentally observed intermediate state of carbon, 
which can be formed from compressing AB stacking 
graphite. Our first-principles calculations based on 
density functional theory show that both the allotropes 
are transparent wide band gap insulators and their 
hardness are comparable to that of diamond. S-carbon 
has the lowest enthalpy compared with bct-C4 carbon, 
M-carbon, W-carbon, and Z-carbon and is the most 
stable new carbon phase theoretically proposed until 
now. H-carbon is always more stable than W-carbon 
and M-carbon. When the pressure is above 5.93 Gpa 
(10.08 Gpa), S-carbon (H-carbon) is more favorable 
than graphite. 



MODELS AND METHODS 

H-carbon and S-carbon belong to the orthorhombic 
lattice and they can be obtained from cold compressing 
AB stacking graphite. Crystal structures of H-carbon 
and S-carbon are shown in Fig. Q] (a) and (d), and their 
corresponding initial graphite structures are shown in 
Fig. [T](b) and (e), respectively. H-carbon with the Pbam 
symmetry containing 16 carbon atoms per crystal cell 
is obtained by reorganizing C atoms to form sp 3 bonds 
through compressing an AB stacking graphite supercell 
with 16 atoms. While, S-carbon with Cmcm symmetry 
containing 24 carbon atoms per crystal cell is obtained 
from larger AB stacking graphite supercell with 24 atoms. 

All calculations are carried out using the density 
functional theory within both local density approxima- 
tion (LDA) [14J . Il5j and general gradient approximation 
(GGA) [16| as implemented in Vienna ab initio simula- 
tion package (VASP) [17|, LLS|- ^ n evaluating of the tran- 
sition pressure from graphite to each superhard phases, 
it is known that LDA is a better choice because LDA can 
give reasonable interlayer distances, mechanical proper- 
ties of graphite sheets due to a delicate error cance- 
lation between exchange and correlation in comparison 
with that of semi-local generalized gradient approxima- 
tion (GGA). The interactions between nucleus and the 
2s 2 2p 2 valence electrons of carbon are described by the 
projector augmented wave (PAW) method [l9|, |20|. A 
plane-wave basis with a cutoff energy of 500 eV is used 
to expand the wave functions. The Brillouin Zone (BZ) 
sample meshes are set to be dense enough to ensure 
the accuracy of our calculations [21]. Crystal lattices 
and atom positions of graphite, diamond, M-carbon, W- 
carbon, Z-carbon, H-carbon and S-carbon are fully opti- 
mized up to the residual force on every atom less than 
0.005 eV/A through the conjugate-gradient algorithm. 
Vibration properties are calculated by using the phonon 
package [22] with the forces calculated from VASP. 



RESULTS AND DISCUSSION 

Similar to M-carbon and W-carbon, H-carbon and 
S-carbon contain distorted five and seven carbon rings, 
as shown in Fig. [T] (c) and (f). Their difference from 
Z-carbon and bct-C4 carbon is the absence of four and 
eight carbon rings. The crystal structure of H-carbon 
belongs to the Pbam space group. At zero pressure, 
the GGA calculated equilibrium lattice constants are 
a=7.874 A, b=4.807A and c=2.524 A. Four inequivalent 
atoms in this crystal occupy the positions at (0.046, 
0.361, 0.500), (0.157, 0.309, 0.000), (0.327, 0.462, 
0.000) and (0.430, 0.393, 0.500), respectively. S-carbon 
belongs to the Cmcm space group and its equilibrium 
lattice constants from GGA calculations are a=2.523 A, 
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FIG. 2: The enthalpy per atom for cubic diamond, hexagonal 
diamond, M-carbon, W-carbon, Z-carbon, H-carbon and S- 
carbon as a function of pressure relative to graphite. 



b=l 1.385 A and c=4.899 A. Four inequivalent atoms 
in S-carbon locate at positions of (0.500, 0.778, 0.250), 
(0.500, 0.442, 0.078), (0.500, 0.132, 0.518) and (0.000, 
0.300, 0.750). Recently proposed Z-carbon contains even 
number carbon rings and hold Cmmm symmetry, with 
equilibrium lattice parameters a=8.772 A, b=4.256 A 
and c=2.514 A. There are only two inequivalent atoms 
in Z-carbon located at (0.089, 0.316, 0.500) and (0.167, 
0.185, 0.000). All the results of the lattice constants of 
diamond, M-carbon, W-carbon, Z-carbon, H-carbon and 
S-carbon derived from both GGA and LDA calculations 
are listed in Tab. HI The LDA calculations tend to give 
smaller lattice constants for all structures considered 
here, but the conclusions derived from LDA and GGA 
are consistent. 

The relative stability of diamond, M-carbon, W- 
carbon, Z-carbon, H-carbon and S-carbon is evaluated 
through comparing their cohesive energy per atom. The 
energy of H-carbon is about 15 meV lower than that of 
W-carbon, about 25 meV lower than that of M-carbon, 
and 10 meV higher than that of Z-carbon at both LDA 
and GGA level. Namely, H-carbon is more stable than 
M-carbon and W-carbon but less stable than Z-carbon. 
Among these allotropes, S-carbon is the most stable 
one whose cohesive energy is about 30 meV lower than 
that of Z-carbon. From the cohesive energy, S-carbon is 
the most stable carbon allotrope theoretically proposed 
until now. The enthalpy per atom for cubic diamond, 
hexagonal diamond, M-carbon, W-carbon, Z-carbon, 
H-carbon as well as S-carbon as a function of pressure 
relative to graphite derived from LDA calculations are 
shown in Fig. [2j The results indicate that H-carbon is 
more stable than graphite when the external pressure 
is larger than 10.08 GPa. H-carbon is always favorable 
than M-carbon and W-carbon at the pressure range 
from to 25 GPa. Amazingly, the transition pressure 
point of S-carbon is only 5.93 Gpa, which is smaller than 



TABLE I: Space group, exchange-correlation functional (EXC), lattice parameters (LP), density (g/cm 3 ), band gap (Eg: eV), 
cohesive energy (Ecoh: eV) and bulk modulus (BO: Gpa) for diamond, M-carbon, W-carbon, Z-carbon, H-carbon and S-carbon. 



Systems Space group EXC 



LP 



Density Eg Ecoh BO 



diamond 


Fd-3m 


LDA 


diamond 


Fd-3m 


GGA 


M-carbon 


C2/m 


LDA a 


M-carbon 


C2/m 


GGA a 


W-carbon 


Pnma 


LDA 


W-carbon 


Pnma 


GGA 


Z-carbon 


Cmmm 


LDA 


Z-carbon 


Cmmm 


GGA 


H-carbon 


Pbam 


LDA 


H-carbon 


Pbam 


GGA 


S-carbon 


Cmcm 


LDA 


S-carbon 


Cmcm 


GGA 



a=b=c=3.536A 
a=b=c=3.574A 
=9.093A, b=2.498A, c=4.108A, beta=96.96° 
=9.194A, b=2.525A, c=4.15lA, beta=97.03° 
a=8.984A, b=4.116A, c=2.498A 
a=9.086l, b=4.156A, c=2.525A 
a=8.677A, b=4.21lA, c=2.489A 
a=8.772A, b=4.256A, c=2.514A 
a=7.792A, b=4.757A, c=2.497A 
a=7.874A, b=4.807A, c=2.524A 
a=2.496A, b=11.293A,c=4.857A 
a=2.523A, b=11.385A,c=4.J 



3.611 
3.491 
3.443 
3.333 
3.455 
3.345 
3.507 
3.394 
3.440 
3.339 
3.489 
3.399 



4.648 
4.635 
3.532 
3.493 
4.325 
4.281 
3.414 
3.273 
4.512 
4.459 
4.451 
4.342 



■8.997 
■7.693 
■8.821 
■7.531 
■8.830 
■7.539 
■8.857 
■7.564 
■8.845 
■7.554 
■8.896 
■7.593 



514.08 
478.96 
447.33 
431.44 
468.12 
447.03 
497.87 
464.13 
466.92 
445.88 
486.29 
468.45 
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FIG. 3: Phonon band structure and density of state for H- 
carbon (a) and S-carbon (b). Electronic band structure and 
density of state for H-carbon (c) and S-carbon (d). 



all of the previously proposed carbon allotropes. To 
further confirm the dynamic stabilities of H-carbon and 
S-carbon, we calculate their phonon band structures and 
phonon densities of state as shown in Fig. [3l Both LDA 
and GGA calculations indicate that there is no negative 
frequency for both structures confirming that H-carbon 
and S-carbon are dynamic stable phases of carbon. 

Space group, density, band gap, cohesive energy 
and bulk modulus of diamond, M-carbon, W-carbon, 
Z-carbon, H-carbon, and S-carbon are also summarized 
in Tab. ID The results indicate that, as super hard 
intermediate phases between graphite and diamond, 
M-carbon, W-carbon, Z-carbon, H-carbon, and S-carbon 



can be formed by cold compressing graphite and their 
densities and bulk modulus are close to that of diamond. 
The density of H-carbon and S-carbon derived from 
GGA is 3.339 g/cm 3 and 3.399 g/cm 3 , respectively, 
which is similar to that of M-carbon (3.333 g/cm 3 ), 
W-carbon (3.345 g/cm 3 ) and Z-carbon (3.394 g/cm 3 ) 
derived from the same method. The bulk modulus of 
H-carbon is 445.88 GPa which is almost the same as 
that of W-carbon (447.03 GPa) and slightly larger than 
that of M-carbon (431.44 GPa). The bulk modulus of 
S-carbon is 468.45 Gpa which is close to that of Z-carbon 
(466.13 Gpa). The values of the bulk modulus indicate 
that H-carbon and S-carbon are superhard materials 
which are comparable to diamond (478.96 Gpa). 

Electronic properties of diamond, M-carbon, W- 
carbon, Z-carbon, H-carbon and S-carbon are also 
investigated under both GGA and LDA calculations. All 
these superhard carbon allotropes are indirect wide band 
gap semiconductors except for S-carbon. The GGA 
results show that H-carbon holds an indirect wide gap of 
4.531 eV which is larger than that of M-carbon (3.523 eV 
), W-carbon (4.325 eV ) and Z-carbon (3.273 eV), and 
smaller than that of diamond (4.635 eV). Different from 
other previously proposed superhard carbon allotropes, 
S-carbon is a direct-band-gap semiconductor with a gap 
of 4.342 eV. The wide band gap indicates that both H- 
carbon and S-carbon are transparent carbon allotropes. 
Fig. 3(c) and (d) show the electronic band structures 
and densities of states of H-carbon and S-carbon derived 
from GGA. The densities of states indicate that the 
valence band maximum for both H-carbon and S-carbon 
is mainly contributed from the 2p states of carbon, 
whereas the conduction band minimum comes from the 
sp hybridized states. 

The simulated X-ray diffraction (XRD) patterns 
for graphite, Z-carbon, H-carbon and S-carbon at the 
pressure of 23.9 Gpa are shown in Fig. |4] to compare 



3 







graphite 

Z-carbon, H-carbon, S-carbon 



13.7 Gpa 



29 (degree) 

FIG. 4: Simulated XRD data for S-carbon, Z-carbon, H- 
carbon, and graphite at corresponding pressure as well as the 
experimental XRD data at 23.9 Gpa and 13.7 GPa for the 
superhard intermediate state derived from reference [5]. 



together with the previous proposed M-carbon, bct- 
C4, W-carbon and Z-carbon will enrich the theoretical 
evidence for understanding the experimental observation. 
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with the experimental data from reference [5]. The 
experimental data under high pressure can be explained 
by H-carbon and S-carbon to some extent. The main 
XRD peaks located in the region between 8.5°- 11° and 
15°- 17° for H-carbon and S-carbon agree well with the 
ex per imental observations. As indicated by Amsler et 
al.[ll| that sole comparison between experimental and 
theoretical XRD can not directly clarify the new carbon 
phase. In view of the cohesive energy, dynamical sta- 
bility, electronic structure, bulk modulus, lowest phase 
transition pressure point, and especially the natural 
transition from AB stacking graphite, H-carbon and 
S-carbon are excellent candidates for the new superhard 
carbon phase. 

Note, added During the course of our submission, 
we became aware of several recent papers [23I430J 
reporting the new superhard carbon phases similar to 
our H-carbon ^ Hfij and S-carbon [11, Hfij 



CONCLUSION 

In summary, we propose two new superhard carbon 
phases, H-carbon and S-carbon, as the possible candi- 
dates for the intermediate phase of cold compressing 
graphite. Both H-carbon and S-carbon with sp 3 carbon 
bonds are optical transparent superhard carbon phases. 
They are more stable than M-carbon and W-carbon, and 
S-carbon is the most table carbon phase theoretically 
proposed until now. Moreover, when the pressure is 
above 5.93 Gpa (10.08 Gpa), S-carbon (H-carbon) is 
more favorable than graphite. These two new members 
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